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ABSTRACT 

We analyze evolution of the basic properties of simulated large scale structure elements 
formed by dark matter (DM LSS) and confront it with the observed evolution of the 
Lyman-a forest. In three high resolution simulations we selected samples of compact 
DM clouds of moderate overdensity. Clouds are selected at redshifts < z < 3 with the 
Minimal Spanning Tree (MST) technique. The main properties of so selected clouds 
are analyzed in 3D space and with the core sampling approach, what allows us to 
compare estimates of the DM LSS evolution obtained with two different techniques 
and to clarify some important aspects of the LSS evolution. In both cases we find that 
regular redshift variations of the mean characteristics of the DM LSS are accompanied 
only by small variations of their PDFs, what indicates the self similar character of 
the DM LSS evolution. The high degree of relaxation of DM particles compressed 
within the LSS is found along the shortest principal axis of clouds. We see that the 
internal structure of selected clouds depends upon the mass resolution and scale of 
perturbations achieved in simulations. It is found that the low mass tail of the PDFs 
of the LSS characteristics depends upon the procedure of clouds selection. 
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1 INTRODUCTION 

At 2 < 1.5 the Large Scale Structure of the Universe (LSS) 
is observed in many galaxy catalogs such as 2dFGRS, SDSS, 
DEEP2, VIRMOS (Percival et al, 2001; Verde et al. 2003; 
Abazajian et al. 2003; Davis et al. 2003; Le Fevre et al. 2005). 
It is manifested as a strong galaxy concentration within rich 
walls and filaments surrounding regions with low density of 
galaxies (voids). At redshifts z > 2 the LSS is observed as 
the Ly-a forest and rare metal systems in high resolution 
spectra of the farthest quasars. The observed properties of 
these populations of the LSS elements are different but now 
it is commonly believed that both of them represent dif- 
ferent manifestations of the same structure formed by DM 
and baryonic components. This statement suggests that the 
Ly-Qf absorbers can be associated with low mass structure 
elements formed by the non luminous baryonic and DM com- 
ponents, while galactic walls and filaments represent fraction 
of the richer LSS elements. This means that the DM struc- 
ture traced by the Ly-o? absorbers is qualitatively similar to 
the rescaled structure observed in the spatial distribution of 



galaxies. In particular, this inference is consistent with ob- 
servations of absorbers in the vicinity of galactic filaments 
and even galaxies within voids (Morris et al. 1993a, b; Pen- 
ton, Stock and Shull 2000, 2002; McLin et al. 2002, Williger 
et al. 2010). 

The facilities of both observational approaches are lim- 
ited but they are complementary to each other. Thus anal- 
ysis of the galaxy surveys is focused on the richer LSS el- 
ements at limited range of redshifts z < 1 — 1.5. It allows 
us to determine the main characteristics of such elements ( 
see, e.g., Doroshkevich et al. 2004) but cannot trace the LSS 
evolution for large redshift intervals. This analysis is concen- 
trated on the investigations of galactic properties and their 
environmental dependence. For example, it allows us to test 
the star formation histories (see, e.g., Panter et al. 2007; 
Skibba et al. 2009) and correlations between orientations of 
the LSS elements and the angular momentum of galaxies 
(see, e.g., Trujillo et al. 2006; Aragon-Calvo et al. 2007; Paz 
et al. 2008; Slozar et al. 2009; Jimenez 2009). 

On the other hand, observations of the Ly-a forest at 
z ~ 2 — 4 allow us to trace some characteristics of the low 
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mass LSS elements and their evolution. However, they rely 
on the observed evolution of neutral hydrogen, which can 
be caused by many factors. Thus, in addition to the actual 
evolution of the DM components of the Ly-a absorbers the 
observed fraction of neutral hydrogen depends upon poorly 
known variations of the UV background. In spite of this, it 
is clear that there is a close link between the observed Ly-a 
absorbers and DM clouds formed at high redshifts. 

It seems that the close link between the Ly-a forest 
and the galaxy LSS can be established by numerical sim- 
ulations of the structure evolution. However, it is not yet 
possible to simulate the LSS evolution in a wide range of 
scales. Thus, simulations performed within large boxes with 
a moderate spatial and mass resolution (for review see Frenk 
2002; Springel, Frcnk& White 2006) reproduce reliably only 
the formation of the richer DM LSS elements similar to those 
observed in the galaxy distribution. On the other hand, high 
resolution simulations performed with the small box sizes 
can reproduce some characteristics of the Ly-a forest (see, 
e.g., Weinberg ct al. 1998; Zhang ct al. 1998; Dave ct al. 
1999; Theuns et al. 1999, 2000; Schayc 2001; Meiksin, Bryan 
and Machacek 2001) but obviously they cannot simulate for- 
mation of the richer LSS elements. Perhaps, more progress 
can be achieved with simulations performed with variable 
resolution (see, e.g., Springel et al. 2008; Diemand et al. 
2008). 

Technical limitations of the box size and resolution lead 
to small and large scale cutoffs of the simulated power spec- 
trum. These cutoffs restrict the potential of such simulations 
and do not allow to reproduce reasonably well the observed 
characteristics of the Ly-a forest (sec, e.g., Meiksin, Bryan 
and Machacek 2001; Gncdin & Hamilton 2002; Tegmark & 
Zaldarriaga 2003; Seljak, McDonald & Makarov 2003; Man- 
ning 2003a,b; Demiaiiski & Doroshkevich 2003). Now such 
simulations are used mostly for the surprisingly stable re- 
construction of the small scale initial power spectrum from 
characteristics of the forest (see, e.g.. Croft et al. 2001, 2002; 
Seljak et al. 2005; McDonald et al. 2005; Viel et al 2004a, 
b). 

An alternative approach is to focus the main attention 
on the regular trends in the evolution of the simulated LSS 
and compare them with observations of the forest properties 
in a wide range of redshifts. Such trends can be revealed 
with special methods applied to representative high reso- 
lution simulations. The successful investigation of the pro- 
cess of halo formation (see, e.g., Navarro, Frenk, & White, 
1995, 1996, 1997; Bullock et al., 2001; Tasitsiomi et al. 2004) 
demonstrates the high potential of this approach. 

The statistical analysis of a large sample of the Ly-a 
absorbers (Demianski, Doroshkevich & Turclianinov 2006; 
Demiaiiski & Doroshkevich 2010) reveals some unexpected 
features of the forest evolution. The most important are 
the regular redshift variations of the mean observed char- 
acteristics of the Ly-a absorbers and the surprisingly weak 
variations of their probability distribution functions (PDFs). 
These results demonstrate the self similar character of evo- 
lution of the observed LSS, which must be tested with nu- 
merical simulations. 

Other important open problems are connected with the 
internal structure of the Ly-a absorbers and its evolution 
with time. They also include an estimate of the degree of 
relaxation of matter accumulated by absorbers with various 



richness, and the link between the thermal and large scale 
bulk motions of the gas and the observed Doppler param- 
eters of the forest. At the present time a reliable observar 
tional discrimination between contributions of these factors 
is problematic. We can also note the surprisingly weak red- 
shift evolution of the observed mean Doppler parameter and 
a complex shape of its PDF indicating the existence of a rich 
sample of DM clouds with small Doppler parameter b < (b). 

Comparison of the observed evolution of the forest with 
the siirmlatod evolution of the DM clouds allows us to clar- 
ify how they are related. For this purpose, in this paper, we 
investigate evolution of the DM LSS in high resolution sim- 
ulations using the Minimal Spanning Tree (MST) and Core 
Sampling approaches (see, e.g., Doroshkevich et al. 2004). 
The space and force resolutions achieved in these simula- 
tions allow us to characterize the evolution of the compact 
DM LSS elements with moderate richness, which is by itself 
an important problem. Thus, we reveal significant differ- 
ences in properties and evolution rates of the LSS elements 
selected in simulations with different mass and space res- 
olutions. These results demonstrate the strong influence of 
small scale perturbations on the internal structure of the 
LSS elements. 

All selected DM elements are significantly larger than 
those observed as the Ly-a forest, what does not allow us 
to perform a direct quantitative comparison of the observed 
and simulated clouds. However, the basic properties of the 
selected sample of DM clouds arc found to be similar to 
the properties of the observed Ly-a forest. In particular, we 
confirm the self similar character of evolution of the DM 
LSS, reproduce the redshift dependence and the PDF of the 
Doppler parameter, estimate the clouds rotation, the contri- 
bution of macroscopic turbulent motions and the degree of 
relaxation of compressed matter measured along the prin- 
cipal axes of clouds. The properties of the LSS elements 
are found to be qualitatively consistent with predictions of 
the Zel'dovich approximation (Zel'dovich 1970; Sliandariri 
& Zel'dovich 1989; Demianski & Doroshkevich 1999, 2004). 

This paper is organized as follows. Simulations that we 
use and methods of their analysis are shortly presented in 
Sec. 2. Evolution of the basic characteristics of the DM LSS 
is described in Sec. 3 and results of the core sampling anal- 
ysis are given in Sec. 4. Discussion of obtained results and 
conclusions can be found in Sec. 5. 



2 NUMERICAL SIMULATIONS 

In this paper we investigate the process of formation and 
evolution of the DM LSS within the concordance cosmolog- 
ical model 

= 0.7, = 0.3, h = 0.7, as = 0.9, n = 1 , 

where Qa and 0„i denote the dimensionless densities of Dark 
Energy and matter, h = Ho/WOkm/s/Mpc the dimension- 
less Hubble parameter, as - the amplitude of density per- 
turbations, n - the power index of the spectrum of pertur- 
bations. 

One simulation - referred further as Siso - was per- 
formed with an MPI version of the Adaptive Refinement 
Tree code (Kravtsov et al. 1997) within a box of Li,ox = 
150/i"^Mpc with 256^ particles, cell size w 0.6/i~^Mpc and 
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mass and force resolutions of 1.7- IO^^/i^^Mq and 18/i~^kpc 
(sec Wojtak et al., 2005 for more details). We have analyzed 
the simulated DM distribution at four redshifts, namely, at 
« = 0, 1, 2 & 3, what roughly covers the interval of observed 
redshifts. 

For comparison we also used the DM distribution ob- 
tained in the MareNostrum Universe (Gottlober, Yepes 
2007). This non-radiative SPH simulation was performed 
with GADGET2 code (Springol 2005) using the concor- 
dance model. This simulation - further referred as 6*500 - 
consists of 1024"^ dark and 1024'^ gas particles in a box of 
Lbox = 500 /i~^Mpc on a side. Both the mass and force 
resolutions (0.75 ■ lO"'/t"^Af0, 15 /i"^ comoving kpc) arc 
comparable to those achieved in the previous sinmlation. 
Such comparison allows us to check the impact of the bary- 
onic component, the code used and the initial realization on 
the properties of the simulated LSS. As a rule, properties of 
clouds selected in this simulation are close to the ones ob- 
tained for the simulation 6*150, what demonstrates the weak 
impact of these factors on the characteristics of the LSS. 

The third simulation - further referred as S50 - con- 
sists of 512^ DM particles in a box of Lbox = 50 h~^Mpc 
on a side. It was performed at the Astrophysical Institute 
Potsdam by M. Steinmetz. With such parameters the ob- 
tained meiss resolution is 7.8 • 10^/i~^Mq and the force reso- 
lution 3/i~^ comoving kpc, which is better than resolutions 
achieved in the previous simulations. However, in this case 
the important influence of large scale perturbations is partly 
suppressed. Comparison of results obtained for these three 
simulations reveals the impact of small scale perturbations 
on the properties of the LSS. In particular, in this Simula^ 
tion we see more complex internal structure of the selected 
LSS elements, earlier formation of high density clouds, and 
we can trace the initial steps of the LSS formation up to 
redshifts « ~ 12 - 14. 



2.1 Selection of high density clouds 

Using the Minimal Spanning Tree (MST) code as described 
in Doroshkevich et al. (2004) we select from the full distri- 
bution of DM particles a set of clouds with three threshold 
parameters that restrict their density and richness. They are 
the threshold overdensity, Sthr, and the minimal and maxi- 
mal richness, Nmin & Nmax, of selected clouds 

<5t/,r = (47r£?hr(Wp>/3)"\ Nrain < Np < Nmax , (1) 

{rip) = (5, 8.6 & 10^)/l^Mpc"^ (p) w 8 • 10^° /i^Mo/Mpc^ . 

Here £thr is the maximal length of the edges of the tree 
within a selected cloud, (p) and (rip) are the mean comoving 
density and number density of DM particles in the simula- 
tions. Using these parameters we can reliably identify clouds 
even of a complex irregular shape. By definition of the MST 
the distances between all neighboring points inside the cloud 
are less than the threshold one, £ < Ithr, and, therefore, 
5 — p/{p) > Sthr- In other words, the selected clouds are 
bound by the surfaces of constant overdensity, 5 = Sthr ■ 

In this paper we are mainly interested in the mildly 
nonlinear evolution of the LSS elements with moderate over- 
density. Such elements dominate at the early period of the 
LSS formation and it can be expected that we observe them 
as the Ly-a forest. Hence the threshold overdensities that 



we adopt in this paper are not very high: Sthr = 1 and 
Sthr = 1-76. 

The other important threshold parameters are the min- 
imal and maximal richness of the selected clouds. The first 
one ensures the reliability and stability of characteristics of 
the selected clouds, while the second one allows us to ex- 
clude from the analysis amorphous multiconnected clouds 
formed at later stages of the LSS evolution. Such extremely 
rich clouds appear owing to the integration of less rich clouds 
and they represent the elements of the well known network 
of the LSS. Evidently, such clouds can not be described by 
just a few local characteristics. 

For these simulations we use the following minimal and 
maximal richness of the clouds under investigation 

iV^in = 30, Nmax = 5000, for Sl50 , 

Nm^n = 70, Nmax = 11 000, for S5OO , (2) 

Nm^r^ = 6.5 ■ 10'', Nmax = lO'' for S50 , 

what approximately corresponds to the same range of cloud 
masses, 

M^i„~ 5- 10" /i-'Mq, M^a^ ~ 8 • 1O"/i-'M0 . (3) 

Further on we refer to the samples selected with 
Mmin < Mcl < Mmax as Samples of compact clouds. For 
some tests we also use samples with Md > Mmin (a full 
sample) and Md > Mmax (a sample of rich clouds). 

The selection of clouds with the core sampling approach 
is described in Sec. 4. 

2.2 Characteristics of clouds 

In this paper clouds are characterized by their principal axes 
determined by their inertia tensor. Such rough description 
allows us to take into account the main global shape of the 
clouds and to estimate their degree of filamentarity and 
sheetness without introduction of any additional parame- 
ters. 

Sometimes much more refined techniques are used for 
selection and discrimination of the DM filaments and walls 
(see, e.g., discussion in Doroshkevich et al. 2004; Aragon- 
Calvo et al. 2007; Zhang et al. 2009). However, advantages 
of such approach are problematic owing to very complex 
shapes, many branches and strong nonhomogeneity of the 
matter distribution typical for both the observed and sim- 
ulated LSS elements. In contrast to our approach such re- 
fined methods are inevitably multiparametric and therefore 
description of so selected clouds becomes more complicated. 
Such approach can be efficient in providing more accurate 
description of local environment of considered objects. 

Here we consider several global characteristics of clouds, 
namely, their mass (richness) and velocity, Md & U, three 
components of their angular momenta Ji and the comoving 
sizes along the three principal axes, Li > L2 > -L3, deter- 
mined through their inertia tensors. 

As is well known (see, e.g, review by Shafer 2008) the 
angular momenta of the clouds are mainly generated by 
the tidal interactions with the surrounding medium (Tidal 
Torque Theory (TTT) Peebles 1969; Doroshkevich 1970; 
White 1984). The angular momentum of a cloud depends 
upon its (random) shape. The mean angular momentum of 
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clouds is small, (Ji) ~ owing to its random orientation 
and (Jf) provides more stable characteristics of clouds. For 
a regular ellipsoidal volume with halfaxes oi > 02 > 03, the 
Zel'dovich theory gives 



( Ji) = 0, (jf ) oc (o^ - alf, i^j^k, 



(4) 



and, so, (J|) > (Ji^s) (Demianski & Doroshkevich 2004). 

In simulations the angular momenta of clouds and the 
contribution of bulk (turbulent) motions are characterized 
by the functions 

Ji = N~'^ ^^tpi, = eifci(Vm - Vc)fc(Xm - Xc)i , (5) 

m=l 

where etki is the totally antisymmetric unit tensor and A'p 
is the number of particles in the cloud. Here Xm & Vm and 
Xc & Vc are the coordinates and velocities of the m*'' particle 
and the center of mass, respectively. 

To characterize the internal dynamics of matter accu- 
mulated by a cloud and its degree of relaxation along the 
principal axes we consider also three velocity dispersions, 
(7i, and three angular momenta of particles, ji, 



(6) 



For a regular ellipsoidal volume with halfaxes ai > 02 > as, 
the Zel'dovich theory gives in contrast to (4) 

(ji) oc {al + aff, i^k^l. (7) 

As is well known, within anisotropic clouds of collision- 
less particles the relaxation occurs independently along each 
principal axis of the cloud. For example, the observed galac- 
tic walls and filaments are almost relaxed along their short- 
est axes but have no time to relax along their longest axis. 
Thus, here we will characterize the degree of relaxation of 
Np particles compressed within a cloud by the parameters 



1 \^ 

Oi = l+ — '^(Vm - Vc)i(Xm - Xc)i/|(Xm - Xc)i| 



(8) 



Here ai is the velocity dispersion and index i shows that we 
consider the relaxation along the i"* principal axis. 

Evidently, the functions Wi discriminate between contri- 
butions of random and regular motions. Indeed, for a cloud 
with relaxed matter distribution and random velocities of 
particles we can expect that jwi — 1| <C 1, because for each 
particle the probabilities to find (v^ — Vc)(xm — Xc)i < and 
(vm — Vc)(xm — Xc)j > are close to each other. In contrast, 
for regular expansion or compression of matter (vr„— Vc) and 
(x,„ — Xc) are correlated and, so, |wi — 1| ~ 1 seems to bo 
much more probable. The case when Wi > 1 or Wi < 1 cor- 
responds to domination of the expansion or compression of 
matter along the i*** axis. 

The core sampling approach was applied only for the 
first simulation and is described in Sec. 4 . 



2.3 Comparison of the simulations with the 
Zel'dovich approximation 

In this paper we consider properties of DM clouds selected 
in the same manner in three simulation. The simulations 



are performed with different codes, in different boxes and 
with different resolution. The cloud selection is performed 
with a complex algorithm and with set of thresholds pre- 
sented in Sec. 2.1. Therefore before we start our analysis we 
estimate the ropresentativity and reliability of the samples 
of selected clouds. For this purpose we compare some in- 
tegral characteristics of selected clouds with predictions of 
the Zel'dovich approximation (Zel'dovich 1970; Shandarin 
&c Zel'dovich 1989; Demianski & Doroshkevich 1999, 2004), 
which reliably describes the mildly nonlinear evolution of 
perturbations. 

For this aim we use here three characteristics of samples 
of clouds, namely, the matter fraction accumulated by the 
LSS, f(M), the PDFs of masses of selected clouds, P(M), 
and the redshift evolution of the velocity of clouds as a 
whole, U(2;). The last test is more sensitive to the large scale 
perturbations, while the first and the second ones are sen- 
sitive to the small scale perturbations, nonlinear processes 
of matter condensation within clouds and the procedure of 
clouds selection. 

The shape of PDFs for clouds sizes and internal veloc- 
ities can be also obtained on the basis of Zel'dovich theory. 
However, here we note only that for these functions and 
the CDM -like power spectrum the Gaussian PDF of ini- 
tial perturbations is transformed to the exponential one (see 
sections 2.4 & 3.3 in Demianski & Doroshkevich (2004) for 
more details). The exponential shape of these PDFs in our 
siirmlations is clearly seen in Figs 1 and 2. In Sec. 3.3 we 
compare predictions of the Zel'dovich theory with the sim- 
ulated angular momentum and turbulent motions (4, 7). 

2. 3. 1 Evolution of the cloud velocity 

The nonlinear processes of matter clustering and relaxation 
of the compressed matter strongly distort velocities of par- 
ticles within clouds but only weakly change the velocities 
of clouds as a whole. As noted in Demianski & Doroshke- 
vich (1999), for the Gaussian initial perturbations, the PDF 
of clouds velocity remains Gaussian with the velocity dis- 
persion closely allied to that described by the linear theory. 
Thus, according to the Zel'dovich approximation the veloc- 
ity dispersion along any principal direction is 



au = 



m = - 



l + z dB{z) 
B{z) dz 



(9) 



1 + 



2.20. 

n 



1 + 1.2Q.r 



-[[l + zf-l], 



= Hl{0^{l + zf + 0^), al 



p{k)dk . 



Here H is the Hubble parameter, Ho = lOOfekm/s/Mpc, 
B{z) describes the growth of perturbations in the linear the- 
ory for ACDM cosmology, p{k) and H{z)as are the power 
spectrum and the amplitude of initial velocity perturbations, 
Q.m and f^A are matter and Dark Energy density parame- 
ters. For both our simulations with the moderate resolution 
we have approximately 
2 07 

as « ^y-^asMpc, B « 1.27[1 + (1 + zf]' 
and for ^; > 1 we get 
H{z) « Ho^M^{l + zf'\ 



-1/3 



B{z) « 1.27(1 -h-?)" 
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au « 353km/s(l + z)'^^'^ {as/0.9) . (10) 

Using the matter velocities in the simulation S150 we 
get for the ID velocity dispersion of clouds at « > 1 

au = 325km/s (1 + ± 0.04) , (11) 

and for the high resolution simulation S50 we get 

au = 315/cm/s (1 + «)"^^^(1 ± 0.04) . (12) 

The differences (~ 8%&11%) between the expected (10) 
and simulated (11, 12) dispersions characterize the impact 
of the large scale perturbations suppressed in simulations 
owing to the finite box size. 

For clouds selected in the simulation Siso with the 
threshold ovordcnsity Sthr = 1-76 and richness S-IO^'^/iMq < 
Mci < 5 ■ 10"/iMq at rodshifts z = 1,2&3 the proba- 
bihty distribution function (PDFs) for the cloud velocities, 
Pu{\Uj\/au), j = 1, 2, 3 are similar to each other with high 
precision and are well fitted by the function 

Pu{xj) = 0.77exp(-xf/2), Xj = \Uj\/au ■ (13) 

Small variations around the Gaussian fit (13) depend on 

pccuharitics of the selected sample of clouds and can also 
indicate the weak mass dependence of the measured disper- 
sions. 

2.3.2 The matter fraction accumulated by the LSS 

To test the simulations and our clouds selection algorithm 
we compare the expected and measured matter fraction ac- 
cumulated by the LSS. According to the Zcl'dovich the- 
ory the fraction of matter expanding along all three ax-es 
is fexp ~ 0.5^ = 0.125 and, therefore, the fraction of matter 
accumulated by all LSS elements is 

ftot - 0.875(1 - z/zf), z^zf. (14) 

This approximate relation slightly overestimates the total 
mass accumulated by pancakes, filaments and clouds. More 
accurate but multipaxametric and more complex estimates 
can be found in Demianski & Doroshkevich (1999, 2004). 
For the simulation S150 we get for Sthr = 1 at ^ < 3 

f{M > 5 • lO"/i"^M0) w 0.73(1 - z/5.7) . (15) 

For the high resolution simulation 5*50 we see earlier matter 
concentration and for Sthr = 1, at z < 6, we get 

f{M > 2.5 • lO^ft^^M©) 0.80(1 - z/10) , (16) 

f{M > 5 • lO"/i"^M0) w 0.67(1 - z/lQ) . 

The moderate differences between the measured (15, 16) 
and expected (14) vahics demonstrate the moderate impact 
of the used selection parameters used, mostly the minimal 
mass. The earlier matter condensation in the high resolu- 
tion simulation S50 (16) as compared with the simulation 
ySiso (15) illustrates the impact of the small scale perturba- 
tions represented by the power spectrum of simulation S50. 

2. 3. 3 The mass function of clouds 

For high density virialized DM clouds the Press-Schechter 
mass distribution is usually reproduced in simulations. How- 
ever, we will see below that for strongly anisotropic partly 



relaxed clouds selected with moderate overdensity the simu- 
lated mass functions are well approximated by a power law. 
Simple qualitative explanation of such distribution is offered 
by the Zel'dovich theory. 

According to this theory the process of the compact 
object formation is determined by three conditions 

Aq,- B{z)AS,^0, i = l,2,3, 

where the function B{z) < 1 was introduced in (9) and 
ASi are the relative displacements of two points separated 
by the comoving distance A^j along i*** principal direction. 
This means that the expected collapsed volume, Vq, and 
corresponding mass, Mzt, can be qualitatively estimated as 

M^t oc Vq = AgiAgaAgs ocVs = AS1AS2AS3 , (17) 

ASi > > 0. 

Here this condition assures that the collapse proceeds along 
all three directions and restricts the volume Vs. 

For any Mo the probability to find the mass of collapsed 
clouds in the range Mo - AM < M^t < Mo + AM is propor- 
tional to the probability to find ASi > Qmin > 0, and the 
volume Vs in corresponding range Vq — AV < 14 < Vb-f AV. 
This simple approach allows us to characterize the shape 
of the PDF P{Vs) oc P{Mzt) but its normalization (in 17) 
requires additional assumptions. More detailed description 
of the complicated process of clouds formation inevitably 
leads to much more complex multiparametric expressions 
(see, e.g., Demiariski & Doroshkevich 2004). 

In turn, the PDF P{Va) can be directly obtained from 
the PDF P{ASi). For this purpose we use the simple three 
step numerical method: 1. to construct the 3D sample of 
ASi in accordance with the theoretical PDF P(ASi); 2. to 
select ID subsample of Vs in accordance with the conditions 
(17); 3. to construct the PDF P{Vs). 

The PDF P{ASi) is Gaussian (Demianski & Doroshke- 
vich 1999, 2004) with dispersion crs(Ag) and weak correlar 
tions, 

iik = {ASiASk)/<T% <!/?,, i^k. 

With this PDF we constructed 10^ realizations of ASi and 
found that the PDF P{Vs) can be fitted by a power law 
P{Vs) oc Vf'. The slope k slightly varies depending on the 
considered range of Vs > Vmin- Thus, for three progressively 
decreasing Vmin = Vi, O.lVi, feO.OlVi the selected subsam- 
ples contain 30 407, 50 025 and 164 103 elements with the 
slop of PDFs K ~ —2.5, —1.8 & — 1.5, respectively. These 
variations indicate that this slop increases for larger Vs . 

In all simulations considered below the PDFs of masses 
of selected clouds are fitted by power laws with k = 
-2, -2.4, -2.2 (19, 20, 22) what is in the range of expected 
values —1.8 > k > —2.5. These variations of k indicate 
mainly the impact of the resolution achieved in simulations 
and the range of selected masses. 



3 PROPERTIES OF SIMULATED DM LSS 

As is well known formation of the LSS begins at high red- 
shifts and its evolution continues up to now. The matter 
fraction accumulated within the LSS elements increases with 
time owing to both the accretion of matter on the earlier 
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formed LSS elements and creation of new low mass LSS el- 
ements and bridges leading to the successive integration of 
the LSS elements into multiconnected richer ones. At small 
redshifts these processes lead to the creation of the joint 
network of compressed matter. The high efficiency of clouds 
integration demonstrates the correlated character of the LSS 
formation and the interaction of small and large scale per- 
turbations. Here we present some quantitative characteris- 
tics of these processes obtained for the LSS formed by the 
DM particles ouly. 

In some important aspects these clouds differ from the 
LSS formed by both the DM and baryonic components and 
observed as the Ly-a forest. In turn, both these structures 
differ from the LSS formed by the luminous matter (galax- 
ies) and observed in large galactic surveys. None the less, all 
these LSS elements are dominated by the DM component 
and their evolution is driven by the same initial perturba- 
tions. So, it can be expected that the most fundamental 
properties of the LSS evolution can be established through 
the analysis of the simulated DM LSS. 

In this section we consider the clouds selected with 
the threshold overdensity, 5thr = 1.76, what allows us to 
partly suppress formation of the very massive multicon- 
nected clouds at redshifts z < 2. Further on we will present 
results obtained for the sample of compact LSS elements 
with DM masses W^^Mq < M < 5 • IO^Mq. 

For the high resolution simulation S50 we can obtain 
more detailed description of the process of matter conden- 
sation within various samples of high density clouds. Thus, 
we have found that within the LSS clouds there is a popula- 
tion of high density cores with 5 > 10. Such early formation 
of high density low mass clouds was found also by Diemand 
et al. (2008). 

3.1 Basic characteristics of compact clouds 

The set of selected compact clouds can be characterized by 
the redshift evolution of their comoving number density, 

(ricisiz)}, masses, {Mci{z)), the fraction of particles accu- 
mulated by clouds, /com, and three comoving sizes {Li{z)) 
{i = 1,2,3) determined along the three principal axes of 
each cloud. 

For the sample of N^a ~ H 000 - 14 000 clouds selected 

in the simulation 5'ir>o with 5thr ~ 1-76 and Mmin < Md < 



Mrr 



have with variations smaller thanlO%: 



(ucUz)) « 3.6- 10 



"^Mpc"^ 



(Mciiz)) 



; 3 ■ 10^^/i"'Mp 



fcom — (j^cls )(Mc,)/(p)«0.13, 
{L,(z))Kl,V^, v = Ma/{Ma(z)), 



(18) 



where (.i are reduced cloud sizes with 

(^i) = 1/r^Mpc, {£2) = 0.5/i"^Mpc, (4) = 0.2h'^Mpc. 

All these parameters characterize the selected sample of the 
LSS elements. Mass dependence of the cloud sizes reflects the 
complex structure of selected clouds and implies their partial 
correlation. Their weak dependence on redshift characterizes 
the balance between creation of new elements and matter 
outflow to the more massive complex LSS elements caused 
by creations of bridges between earlier formed clouds. 

The PDFs of the functions £i plotted in Fig. 1 are also 



IX 0.10 ? 
o : 




Figure 1. The PDFs of clouds sizes £i along the longest, middle, 
and shortest principal axes (top, middle and bottom panels) are 
plotted for redshifts z = 1,2&3 (points, stars and squares) for 
clouds selected in simulation Siso with the threshold overdensity 
5 = 1.76 and richness given by (3). Here Xi = £i/{£i)- 



weakly dependent upon redshift. The exponential decline of 
the PDFs at Xi = £i/{£i) > 1 reflects the rare formation of 
asymmetric objects with large £i while at Xi < 1 the shape 
of the PDFs depends upon conditions used for the cloud 
selection. 

The PDF of cloud masses, Pp{Mci / (Md)) , is surpris- 
ingly independent from redshift but depends upon the range 
of considered velocity dispersion of compressed matter (see 
Sec. 3.5). For all clouds it can be approximated by the power 
law 



Pp{u) = 0.25u, iy = Mci/{Mci) , 



(19) 



what is close to predictions of the Zel'dovich theory (sec. 
2.3.3). It is interesting that here we do not see the exponen- 
tial decline of the PDF, Pp{x), predicted by all theoretical 
models and observed in the luminosity function of galaxies. 
In the simulations such decline appears only at higher red- 
shifts for clouds selected with larger threshold overdensity. 
It can be expected that these differences reflect the com- 
plex structure of richer clouds formed partly by integration 
of earlier formed clouds due to origin of low density bridges 
between them. 

Let us note again that both the mean values (Li) and 
PDFs for the cloud sizes, P{xi = £i/{£i}), plotted in Fig. 1 
do not change with redshift. These results show that for so 
selected clouds, their mean sizes in real space increase with 
time oc (1-1- 2)^^ along each of the three principal axes. This 
means that the mean cloud overdensity does not change and, 
therefore, the possible dissipation of the LSS elements owing 
to the matter expansion is not so essential. 

The characteristics of compact clouds obtained in the 
simulation 5*500 are quite similar to those presented above. 
Let us only note that for this simulation we find some excess 
of low mass clouds as compared with the simulation Siso- 
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Indeed, instead of (19) we get for this simulation the PDF 

Pp{iy) oc V = Ma/ {Ma) , (20) 

what confirms some excess of less massive clouds. This PDF 
differs from (19) but it is in the range of predictions of the 
Zel'dovich theory (sec. 2.3.3) and is quite similar to that 
obtained below for colder less massive clouds (35) . This dif- 
ference shows that some properties of the selected clouds 
are sensitive to such factors as the possible impact of the 
gaseous component, the initial realization of perturbations, 
and the criteria of the clouds selection. 

Analysis of the high resolution simulation S50 allows 
us to obtain more detailed description of the evolution of 
compact clouds. As was noted above in this simulation we 
see formation of the rich population of high density cores 
with the overdensity Sthr > 10 and M^i > 10^° M©. These 
cores have steep density profiles p cx and accumulate up 
to 40% of matter at z = 1. They are usually incorporated 
into less dense clouds. 

For the selected LSS elements we can trace their evolu- 
tion in a wider range of redshifts, 1 < z < 6. Thus, instead 
of (18) we get 

{ricUz)} w 6.3 • 10-^(1 - 0.1;?)Mpc-* , 
(M(z)) w 2.3 • 10^^(1 - O.O7^;);i"^M0, 
fcom ^ 0.18(1 - 0.16;?) , 

L,=£iV^, u = Mal{Ma{z)), (21) 
(^i) = 0.65(1 -I- 0.25a);i"^Mpc , 
(£2) = 0.2(1 -I- 0.5z)h"^Mpc , 
(£3) = 0.1(1 -I- 0.5a);i"^Mpc . 

At 2 = 3 these sizes are quite similar to (18) but at all 
redshifts we see a slow continuous matter compression. This 
difference reflects both the impact of small scale perturba- 
tions and a more detailed description of the LSS evolution 
achieved in the high resolution simulation. 

The PDFs for the cluster sizes and their masses are also 
similar to the previous ones. In particular, 

^» Itl . = M./(M.), (22) 

what is also consistent with predictions of the Zel'dovich 
theory. For larger threshold overdensity the theoretically ex- 
pected exponential decline of the mass function can be seen 
for high z, when majority of clouds are well separated. For 
the mass function of high density halos selected within this 
simulation the exponential decline is clearly seen at all red- 
shifts. 

3.2 The velocity dispersion of compressed matter 

The velocity dispersion of particles within selected clouds 
depends on three factors, which arc: the retained Hubble ex- 
pansion, the enhanced random initial velocities of particles, 
and the velocity created in the course of matter relaxation. 
All these factors are closety linked with the clouds mass 
and therefore the measured velocity dispersion also depends 
upon mass. The relative contributions of these components 
vary with time but for collisionless DM particles the rate of 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 



X 

Figure 2. For the simulation 5150 the PDFs of the reduced ve- 
locity dispersion, ci_2.3, along the longest, middle and shortest 
principal axes (top, middle and bottom panels) are plotted for 
redshifts z = 1,2&;3 (points, stars and squares) for clouds se- 
lected with the threshold overdensity 5 = 1.76 and richness given 
by (3). Here Xk = Ck/{ck), A; = 1,2, 3 . 

these variations depends upon degree of matter compression 
and relaxation achieved along each of the principal axes. 
This means that the velocity dispersions along these axes 
are different and must be determined separately along the 
longest, oi, middle, CT2, and shortest, CT3, principal axes of 
each cloud. 

It is important that these velocity dispersions are much 
smaller than the dispersions au introduced in (11). This 
difference is caused by the strong correlation of the parti- 
cle velocities at the scales of clouds. The simple analysis 
shows that the dispersion of initial velocities within clouds 
decreases with the cloud mass, what stimulates the formar 
tion of low mass high density clouds. 

For clouds selected in the simulation 5i5o with 5thr = 
1.76 and Mmin < Md < Mmax the velocity dispersions 
along three principal axes, at, slowly vary with time but 
they depend upon the cloud richness, 

ai = CiZ/"% ai = 0.5, a2 = 0.4, aa = 0.3 . (23) 

The mean values of the reduced velocity dispersions Cj are 

(ci) = 58/cm/s, (C2) = 66A;m/s, (cs) = 80km/s , (24) 

and their PDFs, P{xi = Ci/(ci)), are plotted in Fig. 2 for 
three redshifts z = 1, 2, 3. The very weak redshift evolu- 
tion of these PDFs verifies the self similar character of the 
simulated LSS evolution. 

These PDFs demonstrate the expected exponential de- 
cline P{xi) oc exp(— a;i/0.3) for Xi > 1. At the same time, 
there are cold clouds with a < (d). 

For the simulation S^oo the velocity dispersions, (c;) are 
~ 20% smaller than (23). It is consistent with the excess of 
low mass clouds noted above (20) and confirms the sensi- 
tivity of the characteristics of the low mass clouds to the 
mass resolution, impact of the gaseous component and real- 
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Figure 3. The PDFs of the reduced integral angular momentum, 
/i, and the reduced integral intensity of turbulent motions, T, are 
plotted for redshifts z = 1,2&3 (points, stars and squares) for 
clusters selected with the threshold overdensity 5 = 1.76 and 
richness given by (3). Here xi = n/ {n), X2 = t/{t). 



ization of the initial velocity field. The PDFs of the velocity 
dispersions are also similar to those plotted in Fig. 2. 

For the simulation 550, for the same mass range, at red- 
shifts 1 < z < 5 the mean velocity dispersions are described 
by the relation (23) with 

ai = 0.4, q;2,3 — 0.3 

(ci) ^ 105(1 - Q.lbz)kni/s, (25) 
(C2) ~ 95km/s, (cs) ^ SOkm/s. 

Their PDFs are quite similar to those plotted in Fig. 2. 

Differences between estimates (23, 24) and (25) can be 
related to the influence of small scale perturbations result- 
ing in a more rapid relaxation and more complex internal 
structure of selected compact clouds. This is seen as the for- 
mation of high density cores within these clouds and as a 
decrease of power indices ai &a2,3 in (25) as compared with 
(23). In turn, action of these factors increases the reduced 
velocity dispersions C1&C2 (25) as compared with (24). At 
the same time, for both simulations 03 and C3 remain quite 
similar. 



3.3 Rotation and turbulent motions of the 

compressed matter 

In the partly relaxed clouds particles move on closed orbits 
what generates both the turbulent motions, ji (6) and the 
angular momenta of clouds as a whole, Ji (5) along three 
principal axes. Evidently, for isotropic simulations (Ji) = 
and further on we will consider characteristics of the abso- 
lute values of angular momenta \Ji\ only. As before, these 
functions depend upon cloud richness. For clouds selected in 
the simulation we have 



\Ji\ 



(26) 



where /ii & Ti characterize dimcnsionlcss reduced angular 
momenta and intensity of the turbulent motions, Li & ai 
(18&23) are the size and velocity dispersion of clouds. The 
mean values of the reduced angular momenta are found to 
be independent of z and 



(mi) = 0.08, (m2) = 0.15, (H3) = 0.15 . 



(27) 



The small clouds rotation around the longest axis and ex- 
ponential shapes of PDFs 



P(a;i) w 0.12exp(-a;i), Xi=m/{m}, 



(28) 



are consistent with the TTT predictions. 

The mean values of the reduced intensity of turbulent 
motions are also found to be independent of z and 

(n) = 0.64, (r2) = 0.76, (ra) = 0.65 , (29) 

and their PDFs are fitted by the functions 

P(a;i) «0.14exp[-(a;i-l)V0.14], x^ = Ti/{Ti). (30) 

For the total momentum of the clouds, J, and for the 
total intensity of turbulent motions, j, 



gl = (L2 + -^'3)(o-2 + o-f). 



92,3 = (Ll + L3 2){Cri + a3^2) 



J = \/Jl + H + H = Mt, 3 = Vjf+ji+js = Tgt , (31) 
Qt = (L? + + Ll){af + a2+ 0-3) , 
we get 

(m>=0.17, (r) = 0.8. (32) 

The PDFs for the functions /i&r arc plotted in Fig. 3. The 
obtained value of (r) corresponds to domination of elliptical 
trajectories of particles with the ratio of axes b/a ~ 0.57, 
what is consistent with the ratios of the cloud sizes (18). 

For the simulation S50 we have, at 1 < z < 6 and for 
the same mass range, instead of (27), (29) and (32): 

im) = 0.08(1 - 0.07«), (n) = 0.57, 

(ai2) = 0.1(1 - 0.07«), {t2} = 0.60(1 - O.OSz), (33) 

ins) = 0.09(1 - O.lz), {t3} = 0.57(1 - 0.09^) , 

in) = 0.13(1 - 0.08«), (r) = 0.7(1 - 0.08^) . 

Comparison of (fi) with (r) in (27), (29), (32) and (33) 
shows that only ~ 15 — 20% of the angular momenta of indi- 
vidual particles are transformed into the angular momentum 
of a cloud as a whole, what characterizes the efficiency of the 
generation of angular momentum according to the TTT. 

These results are also consistent with predictions of the 
Zel'dovich theory (4 - 7) and indicate the significant contri- 
bution of initial random velocities in the turbulent motions 
within compressed clouds. 

3.4 Relaxation of the compressed matter 

The three parameters wi,W2&i ws introduced in (8) allow us 
to discriminate between the regular expansion and compres- 
sion of matter along the principal axes and to estimate the 
degree of matter relaxation. As noted in Sec. 2.2, > 1 
indicates domination of matter expansion, while < 1 
corresponds to domination of matter compression along the 
i*'' principal axis. For the simulation Sibo and for redshifts 
z = 1,2,3 the mean values (wj) are listed in Table 1 and 
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Table 1. Degree of relcixation of selected clouds 





z=l 


z=2 


z=3 


(wi)-l 


0.05 


-0.2 


-0.42 


{W2)-1 


-0.4 


-0.53 


-0.62 


{w:,)-l 


-0.2 


-0.42 


-0.55 


frel 


0.8 


0.44 


0.2 



their PDFs are plotted in Fig. 4. These results indicate that 
for our samples of selected clouds the matter compression 
along all three axes dominates but it is progressively decel- 
erated with time. The symmetry of the PDF P{wi) points 
out that for our selected clouds both the compression and 
expansion along the longest axis are almost equally typical. 
These estimates only weakly depend upon the richness of 
clouds. 

Another manifestation of the retained expansion along 
the longest axis is the tight correlation of velocity dispersion 
and size: 

rO.8 

with the deviation of ~10%. This expression is quite similar 
to the Hubble flow. For the middle and the shortest axes 
such correlations arc very weak. 

These results confirm that along the longest principal 
axis the relatively slow rela:xation is accompanied by the re- 
tained expansion of clouds, which persists for some time 
after cloud formation. In contrast, along the middle and 
shortest principal axes the moderate retained compression of 
clouds predominates. These results are consistent with the 
strongest matter compression along the shortest principal 
axis (which is responsible for pancakes formation) and the 
successive transformation of prolate pancake - like clouds 
into elongated filamentary - like ones. 

To estimate the degree of matter relaxation along the 
shorter axis we can use the fraction of clouds, frei{z), for 
which |1 — W3(2;)| < 0.5. For three redshifts the estimates of 
this fraction arc listed in Table 1. The rapid growth of this 
fraction with time illustrates the progressive relaxation of 
the compressed matter along the shortest principal axis. The 
fraction of matter relaxed along the longest principal axis is 
only half of that along the other two axes. For the simulation 
Ssoo the functions {wi{z)) — 1 and fr{z) are similar to the 
ones discussed above. 

More complex internal structure of clouds selected from 
the high resolution simulation Sso leads to a more rapid 
relaxation of the compressed matter. Thus, we have for the 
matter fraction with 11 — wsl < 0.5 



frel (.? = 5) « 0.6, freliz = 3) ^ Q.9 . 



(34) 



These results illustrate the impact of small scale perturba- 
tions and strong links between the resolution achieved in 
the simulation and the processes of clouds formation and 
relaxation. 



3.5 Colder and hotter clouds 

The strong link between the velocity dispersions and rich- 
ness of clouds (23) illustrates the large scale correlation of 
the initial velocity field. None the less, the scatter around 
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Figure 4. For the simulation S150 the PDFs of the ratio of pro- 
jected velocity and velocity dispersion along tiic longest, wi , mid- 
dle, W2, and shortest 103, axes (top, middle and bottom panels) 
are plotted for redshifts z = 1,2&3 (points, stars and squares) 
for clouds selected with the threshold overdensity <5 = 1.76 and 
richness given by (3) 



the mean dispersions (23) is large and differences between 
properties of colder and hotter clouds deserve special inves- 
tigation. 

In the simulation 5*150 the subpopulation of colder 
clouds was selected from the full sample of clouds by the 
condition that as < 40km/ s. It contains about half of 
clouds with the lower richness and the steeper mass func- 
tion Pc(Me(/(Me()): 



10" ft" ^ Mb, 



Pc{Mcl) OC [Mc(/(Me()]- 



(35) 



The subpopulation of hot rich clouds with aa > 
40km/ s, is characterized by larger mean richness and shal- 
lower mass function, Ph{Mci{z) / {Mci{z)}) , 

{M,i{z)) « 3.2 • lO^^ft-^Mo, Ph{M,i) OC [M,i/{M,i)\-^-^ . 



(36) 



Both colder and hotter subpopulations are dominated 
by partly relaxed clouds expanded along the longest prin- 
cipal axis and compressed along the middle and shortest 
axes and for them the numerical estimates of (wi) are quite 
similar to the ones listed in Table 1. For the subpopula- 
tions of hotter clouds the fraction of relaxed clouds with 
|1 — ^3(2;)! < 0.5 is close to frei listed in Table 1, while for 
the colder clouds this fraction is smaller by a half. 



4 CORE-SAMPLING APPROACH. 

For a more detailed comparison of the simulated matter dis- 
tribution with the observed characteristics of the LSS and 
the Ly-a forest it is convenient to use the core-sampling 
approach (Doroshkevich et al. 2001, 2004; Dcmiahski ct al. 
2006). With this approach we retain only particles within 
the selected sample of clouds, and divide the simulated 
box into a system of rectangular cores of comoving size 
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dcore X dcore X L^ox, and considei properties of particles sit- 
uated within each core. 

All particles within a core are projected on the core axis 
and using the ID cluster analysis are collected into clumps. 
For so selected clumps wc can determine several characteris- 
tics such as the mass, Mcr, and the comoving surface density, 
qcr = Mcr/d^ore- Other characteristics relate to the clump 
position and velocity along the core. They are the comov- 
ing clump separation, Dsep, and their size (thickness), tcr, 
and overdcnsity, Scr = Qcr/tcr/{p), the relative velocity of 
neighbouring clumps, \5U\, the dispersions of clump veloc- 
ity, ffu, and particle velocity within clumps, cr„. We can also 
estimate the degree of matter relaxation, Wc,-, (8), and the 
dispersion of transverse velocity of matter within clumps, 

The core-sampling analysis utilizes four parameters, 
namely, the comoving size of the core, dcore, the comov- 
ing threshold linking length used in the ID cluster analy- 
sis, Ithr, and the minimal and maximal richness of clumps, 
qthr & Qmax, retained for the analysis. Results of the analysis 
depend mainly upon the choice of dcore- Two last parame- 
ters allow us to remove from the sample poor clumps that 
arise owing to random intersection of the core and a cloud 
periphery, and small number of extremely rich clumps with 
q > Qmax, what strongly distorts the mean measured char- 
acteristics. 

Here we apply the core-sampling approach to investi- 
gate the full sample of clouds selected in the simulation Siso 
with the threshold overdensity Sthr = 1.76. In our analysis, 
we use 

dcore = l.hr^MpC, Ithr = 0.5/i"^Mpc, (37) 
qthr = Mcs/diore = 3 • lO^^Mo/lMpc"^, Qmax = Wqthr ■ 

These samples contain ~ (1.3— 1.5) • 10^ clumps. Naturally, 
variations of these parameters change the fraction of low and 
high mass clumps in the selected sample and also the main 
characteristics of the sample of clumps. However for a wide 
set of limits (37) we see a weak redshift dependence of the 
mean characteristics and their PDFs, what confirms the self 
similar character of the LSS evolution. 

In many respect this approach and its limitations are 
similar to those appearing in the analysis of the Ly a forest 
(Demianski et al. 2006) and the pencil beam observations 
of the LSS (Doroshkevich et al. 2001, 2004). Of course the 
criteria of clump selection in observations and simulations 
are not identical. Indeed, if the Doppler parameter of the 
forest line coincides with the velocity dispersion a„, then 
the evolutions of the DM surface density q{z) and the col- 
umn density of neutral hydrogen, Nhi, are very different as 
they are driven by different factors. Thus, we cannot directly 
compare simulated results with observations. However even 
such analysis allows us to reveal qualitative similarity of the 
observed and simulated LSS evolution and to demonstrate 
the impact of the main factors that influence evolution of 
the observed LSS. 

4.1 Dynamical characteristics of the DM LSS 

The redshift evolution of dynamical characteristics of the 
simulated DM clumps obtained with the core sampling ap- 
proach are found to be similar to those obtained in the pre- 
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Figure 5. The PDFs of clusters velocity, P(xu), the velocity 
dispersion within clusters, P(xv), and the degree of relaxation, 
P(wcr); arc plotted for redshifts z = 1,2&3 (points, stars and 
squares) for clusters selected within cores under conditions (37). 
Here Xu = \U\/au, Xy = <t„/((Tv). 

vious section. Thus, with the threshold parameters (37) and 
for 3 > > 1 we get for the velocity dispersion of clumps, 
cT(7, and for the relative velocity of neighboring clumps, 5U 

{au} w 308A;m/s(l + 2)"^/^(l. ± 0.05) , (38) 

{\6U\) w 272/cm/s(l + x)"^''^(l. ± 0.03) . (39) 

So defined au and the value (11) measured for all particles 
are close to each other. The PDF for the dispersion of the 
ID clumps velocity, P{xu = \U\/au), plotted in Fig. 5 (top 
panel) is close to the Gaussian one. However at all redshifts 
there is some excess of clumps (~ 6%) with extremely low 

(Tu- 

Both velocity dispersions of the compressed matter 
measured with the core sampling approach along the core, 
((T„), and in the transverse directions, (<Ttr), are weakly de- 
pendent upon redshift, what is similar to the behavior of the 
cTj (23). Thus we get 

{ffv} w {crtr} w 49A;m/s(l ± 0.2) . (40) 

The random orientation of the cloud principal axes and 
the probing core leads to intermixture of velocities along 
three principal axes, and the velocity dispersion measured 
in transverse directions, {fftr}, is the same as {cr„) measured 
along the core. 

The PDFs of the velocity dispersion of matter com- 
pressed within clumps, P{xv ~ o-i,/{(t„)), is plotted in Fig. 
5 (middle panel) for three redshifts, z — 1,2, 3. As is seen 
from this Figure for Xv > 0.5 these PDFs arc quite similar 
to each other and they are well fitted by the exponential 
function 

P{xv) w 1.5ea;p(-a;„/0.75) . (41) 

However fraction of clumps with a;„ < 0.1 progressively de- 
creases with time from ~ 25% at z=3 and down to 10% 
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at redshifts z = 1. This fraction can be related to the low 
mass unrelaxed clumps formed due to the intersection of 
cores with complex low density periphery of the selected 3D 
clusters. 

As before, the mean value (wcr) ~ 1 characterizes the 
symmetry of the matter compression and expansion along 
the core in the sample of selected clumps. However, the PDF 
P{wcr) plotted in Fig. 5 (bottom panel) characterizing the 
degree of matter relaxation along the core shows some vari- 
ation with redshift. Using this PDF we can estimate the 
fraction of clumps, frei{z) with 0.5 < Wcr < 1.5 which are 
almost relaxed along the core at redshifts 2: = 3, 2, & 1 as 
follows: 

freli3) = 0.56, /,ei(2) = 0.64, frcl(l) = 0.7. (42) 

These fractions are similar to the fraction of clusters relaxed 
along the shortest axis listed in Table 1, what demonstrates 
moderate impact of the macroscopic motions along the core. 
This result seems to be natural because of the moderate 
difference between the fraction of matter relaxed along the 
shortest and longest principal axes (Table 1& Fig. 4). 

4.2 Spatial characteristics of the DM LSS 

The simulated overdensity of clumps, (Scr), the surface den- 
sity, {?cr(l + z)^^), and the comoving thickness of clusters, 
(tcr), are found to be weakly dependent upon the redshift, 

(Ml + ^n«4.10"^(l.±0.2), 

(fcr) « 1 /i"^Mpc(l. ±0.1) , (43) 

where (p) is the mean comoving density (1). 

The redshift evolution of the comoving cluster separa- 
tion, 

{D,ep) » 16.2h-'Mpc(l. ± 0.1) , (44) 

is also weakly dependent upon redshift but it strongly 
depends upon threshold parameters (37) determining the 
clumps selection procedure. 

For the simulation S50 the mean comoving separation 

between clumps is 

(D.ep) ~ 7/r'Mpc. (45) 

The difference between this value and (44) illustrates the 
impact of the box size (the limited core length). The same 
result as (45) is obtained when we cut a "box" of 50ft~^Mpc 
from the simulation 8150- 

The PDFs of the DM surface density, P{qcr/{<lcr)), and 
of the objects separation, P{Dsep/{Dsep)), also weakly de- 
pend on redshift and are well fitted by exponential functions 

P{xq) w 0.8exp(-a;q/0.8), Xg = qcr/{qcr) > 0.2 , (46) 

P(^Xsep^ ^ 0.9 exp( 3'sep); Xsep — P^sepj {^P^ sep) • 

These results show that with the core sampling ap- 
proach the redshift evolution of the simulated LSS is well 
characterized by the slow regular variations of its mean 
characteristics, while the corresponding PDFs practically 



do not change at all. This fact verifies again the self sim- 
ilar character of the LSS evolution, at least for its richer 
elements formed by the DM component of the Universe. It 
is consistent with the observed evolution of the Ly-a forest 
(Demianski et al. 2006). 



5 SUMMARY AND DISCUSSION 

Many branches and strongly nonhomogeneous matter distri- 
bution typical for the richer LSS elements determine their 
complex multiconnected structure and do not allow to char- 
acterize them in any simple way. Thus, we are compelled to 
restrict our analysis to the sample of compact DM clouds 
with a moderate richness and overdensity. The simple pro- 
cedure of the sample selection and the basic characteristics 
of clouds are described in Sec. 2. More refined technique pro- 
posed for selection and discrimination of the DM filaments 
and walls (see, e.g., discussion in Doroshkevich et al. 2004, 
Aragon-Calvo et al. 2007; Zhang et al. 2009) can be efficient 
mainly for the solution of some special problems. 

In many respects evolution of the selected population is 
typical for all LSS elements but some discussed characteris- 
tics depend upon the procedure of sample selection. First of 
all this relates to the spatial characteristics of clouds. 

As is well known the LSS evolution is determined by 
the accretion of diffuse matter, creation of poorer clouds 
and integration of clouds into larger ones. These processes 
lead to rapid matter concentration within population of 
large multiconnected clouds and to successive network for- 
mation. Thus, in the simulation 5*150 the matter fraction 
accumulated by richer clouds with M^i > lQ^*'h^^ Mq in- 
creases by a factor of 5 between redshifts = 3 and z = Q, 
while the fraction of particles accumulated by clouds with 
5 • lO"/i"^M0 < Np < lO"/i"^M0 increases by a factor 
of 1.5 onlj'. This means that at z < 3 the LSS evolution 
is dominated by successive integration of the earlier formed 
LSS elements into the richest ones. The evolution of compact 
clouds with moderate richness is determined by a balanced 
action of coalescence and formation of new clouds. 

The existence of the population of DM clouds with mod- 
erate richness is in itself quite interesting. In many respects 
its slow evolution leads to the stability of measured clouds 
characteristics, which are distorted mainly by the successive 
formation of high density cores. However, at higher redshifts 
the infiuence of the process of integration of clouds becomes 
less significant, the creation of new clouds and accretion of 
diffuse matter dominate and characteristics of clouds (21) 
evolve more rapidly. 

The performed analysis of the simulated DM LSS al- 
lows us to clarify some important factors that cause its evo- 
lution, and to link them with the power spectrum of the 
initial perturbations and with the observed evolution of the 
Ly-a forest. The limited mass resolution - M ~ 2 ■ lO^^AfQ 
and M ~ 7 ■ W"^ Mq in analyzed simulations - prevents the 
formation of low mass clouds comparable with the majority 
of the observed Ly-a absorbers. This means that our results 
cannot be directly compared with the observed characteris- 
tics of the Ly-a forest. However, such analysis allows us to 
find some characteristics of evolution of the simulated DM 
LSS elements, what is very interesting in itself. It is impor- 
tant also that the main basic features of this evolution are 
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similar in many respects to tiie observed properties of the 
Ly-a forest, what allows us to explain some peculiarities of 
the forest evolution discussed in the Introduction. 

5.1 Main inferences 

The most important features of the evolution of the selected 
DM LSS elements can be summarized as follows: 

(i) Weak redshift variations of the basic characteristics 
of the selected sample of DM clouds such as their mean 
comoving size and velocity dispersion. 

(ii) Measured PDFs of the basic characteristics of the LSS 
elements are weakly dependent upon the redshift, what im- 
plies the self similar character of the LSS evolution. 

(iii) Significant degree of relaxation of the compressed 
matter along the shortest principal axis and retained matter 
expansion along the longest principal axis. 

(iv) Weaker rotation of clouds along their longest axis. 

(v) Some measured clouds characteristics depend on the 
procedure of clouds selection. 

5.2 Impact of simulation parameters 

Comparison of characteristics of clouds obtained for three 
simulations with different box sizes and resolutions shows 
that some of them are sensitive to the resolution, the box 
size and the impact of baryonic component. In particular, 
complex internal structure of clouds selected in the high res- 
olution simulation S50 leads to more complex evolution of 
the basic characteristics of clouds as compared with the re- 
sults obtained for lower resolution simulations. In the simu- 
lation Sr,o clouds arc partially fragmented into high density 
cores, what distorts their dynamical properties. The fraction 
of high-density virialized clumps increases with resolution 
achieved in simulation. However, even a large mass fraction 
(up to 40%) of these high-density clumps does not affect 
strongly the main qualitative inferences: the self-similar evo- 
lution of characteristics of clouds and exponential decline of 
their PDFs. 

The most prominent impact of resolution is detected on 
the clouds size along the middle and the shortest axes £2,3 
(21), the velocity dispersion along the longest principal axis 
(25), and the relaxation parameter, while other properties 
of clouds are quite similar in all three simulations. 

In the simulation 5*50 the box size affects separation be- 
tween clouds (44, 45). The box size of 150 Mpc should 
be enough for the stable reproduction of both these proper- 
ties. 

5.3 Self similarity of the DM LSS evolution 

The most interesting result of our analysis is the self similar 
character of evolution of the DM LSS which is manifested 
at 2 < 4 — 5 as the surprising stability of PDFs of the basic 
characteristics of clouds found with both the 3D analysis 
in Sec. 3 and with the core-sampling approach in Sec. 4. 
This result is obtained for three simulations performed in 
different way and with different rmmber of particles, reso- 
lution and realization of the initial particle positions and 
velocities. These features are similar also to the redshift in- 
dependence of PDFs of the observed characteristics of the 



Ly-a forest (Demianski et al. 2006). These analogies point 
to the universal character of the detected self similarity. 

Other manifestation of this self similarity is the uni- 
versality of the NFW density profiles within simulated DM 
halos (see, e.g., Navarro, Frcnk, White 1995, 1996, 1997). 
The self similar character of the DM halos formation is seen 
as the regular redshift variations of the internal structure of 
individual halos. However, the complex process of formation 
and relaxation of the high density LSS elements and halos 
can bo described analytically as a self similar process only 
for the simplest cases (see, e.g., Fillmore & Goldreich, 1984; 
Gurevich, & Zybin, 1995; Sikivie, Tkachev, Wang Yun, 1997; 
Nusser, 2001). 

For objects formed in the course of mildly nonlinear 
compression the self similar character of evolution is pre- 
dicted by the Zel'dovich approximation (Zel'dovich 1970; 
Shandarin & Zel'dovich 1989), where the evolution is de- 
scribed by a product of time and space dependent func- 
tions. The statistical description of the DM LSS based on the 
Zel'dovich theory (Demianski & Doroshkevich 1999, 2004) 
clearly demonstrates the expected self similarity. 

The simulated evolution of DM clouds is a strongly de- 
terministic process driven mostly by the used realization of 
initial perturbations. The statistical properties of clouds are 
determined by the correlation and structure functions ex- 
pressed through the simulated power spectrum. Therefore 
we can expect that the self similarity of the simulated evo- 
lution is caused by the properties of the power spectrum 
responsible, in particular, for the formation of clouds. 

Such self similarity naturally appears for the power 
spectrum approximated by a power law. Indeed, in this case 
all correlation and structure functions are described by the 
time dependent amplitudes and standard space dependent 
functions. This usually leads to the self similar evolution 
of the basic statistical characteristics of the DM LSS. The 
power spectra used in considered here simulations arc well 
approximated by a power law, what explains the found self 
similarity. 

Of course, this inference is related to the limited range 
of scales of the power spectrum and therefore to the limited 
range of clouds richness. However, we can expect that this 
range extends down to the scales typical for the majority of 
the forest elements, what in turn explains weak redshift vari- 
ations of the observed PDFs of forest characteristics. In con- 
trast, the more complex process of formation of the largest 
clouds and the network of the high density LSS elements is 
driven by the large scale part of the power spectrum that de- 
termines the formation of low density bridges between high 
density clouds. It is convenient to analyze these mildly non- 
linear processes on the basis of the Zel'dovich approximation 
(Demianski & Doroshkevich 1999, 2004). 

5.4 Evolution of the velocity dispersions and the 
degree of relELxation of clouds 

Second important result of our analysis relates to the evolu- 
tion of the internal structure of clouds. Both the moderate 
mass resolution achieved in simulations and moderate rich- 
ness of the dominant population of clouds prevents more 
detailed description of the anisotropic matter distribution 
within clouds. However, high force resolution achieved in 
the considered here simulations allows us to estimate, with 
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a reasonable precision, tlie evolution of the velocity disper- 
sion and the degree of rela^cation of the compressed matter. 

The mean values of the velocity dispersions cri , (T2 & era, 
(23) and even their PDFs presented in Fig. 2 for three red- 
shifts arc quite similar to each other. However, the PDFs 
of the degree of rela^Kation presented in Fig. 4 demonstrate 
significant differences in characteristics of the velocity com- 
ponents. Indeed, at all redshifts the excess of larger values 
of wi > 1 indicates noticeable contribution of the retained 
matter expansion and moderate degree of relaxation along 
the longest ax;is. At the same time, the excess of smaller val- 
ues of W2 < 1 demonstrates domination of matter compres- 
sion along the middle axis and progressive transformation of 
a pancake like cloud into a filamentary like one. The degree 
of relaxation along this axis is moderate but it is larger than 
that along the longest axis. 

In contrast, the more symmetric shape of the PDF 
P{w3) in Fig. 4 indicates that similar fraction of matter 
is expanded and compressed along the shortest axis. This 
Figure shows progressive point concentration near the cen- 
ter, what refiects the progressive randomization of velocities 
and growth of the degree of relaxation of the compressed 
matter. Evaluation of the fraction of relaxed clouds and the 
corresponding matter fraction frei (Table 1) characterizes 
these processes quantitatively. 

Similar characteristics of the velocity dispersions are 
also obtained with the simulations S500 and 850- The most 
serious quantitative divergences arc found for characteristics 
of low mass clouds and are discussed in Sec. 3. 

It is interesting, that the estimated state of relaxation 
of clouds with the core sampling approach in (42) is very 
close to that obtained with the 3D analysis. This fact shows 
that the random orientation of the clouds with respect to 
the line of sight results in the intermixture of thermal and 
large scale bulk velocities but does not distort strongly the 
measured velocity dispersions and the degree of relaxation 
of compressed matter. 

Special analysis confirms that the subpopulations of 
colder and hotter clouds with as < (as) and as > (as) re- 
spectively, contain comparable fractions of matter with dif- 
ferent mean richness and mass functions. Both subpopula- 
tions are dominated by partly relaxed clouds slowly expand- 
ing along the longest principal axis. The colder clouds rep- 
resent population of anisotropic low mass clouds and both 
the shape of their mass function and PDFs of other char- 
acteristics arc strongly depended upon parameters used for 
selection of the considered DM clouds. In contrast, for popu- 
lation of richer clouds the PDFs of cloud characteristics only 
weakly depend on the parameters that select these clouds. 

A strong difference between velocity dispersions au (9) 
and ai (23) related to the motion of clouds as a whole and 
to velocities within clouds is caused by the large scale spa- 
tial modulation of the initial velocity field in regions of the 
clouds formation. The same effect is seen in estimates of the 
relative velocities of neighboring clusters, AU (39). 

5.5 Angular momentum of the LSS elements. 

Recently the angular momenta of the LSS elements arc 
widely discussed and compared with the alignments and an- 
gular momenta of galaxies and their DM halos (see, e.g., 
Vitvitska et al. 2002; Patiri et al. 2006; Lee & Erdogdu 2007; 



Park et al. 2007; Aragon-Calvo et al. 2007, 2008; Cucsta 
et al. 2008; Paz et al. 2008; Slosar et al. 2009; Jimcncs ct 
al. 2009; Zhang et al. 2009; Lovell et al. 2010). The main 
inferences of the TTT and, in particular, preferential orien- 
tation of the angular momenta along the shortest axes of 
the LSS elements are basically confirmed. Infiuence of the 
anisotropic matter infall into the LSS elements on their an- 
gular momenta was discussed by Doroshkevich (1973) for 
baryonic component and in Bernardo et al. (2002) for DM 
particles. However, analysis of simulations shows that for 
DM component contribution of this effect is moderate and 
our results obtained in Sec. 3.4 are consistent with TTT 
predictions. 

On the other hand, large scatter is found between the 
shape of the host LSS element and the angular momenta of 
galaxies and their DM halos incorporated in that element. 
This scatter is expected and can be caused by action of 
several factors. 

Thus, comparison of the functions J&ij (31) shows 
that the process of generation of angular momentum of DM 
clouds is only moderately efficient (less than 20%). This fact 
in itself suggests a weak correlation between the angular mo- 
menta of the host LSS clement and high density low mass 
subclouds and it basically explains the scatter mentioned 
above. Moreover, more detailed analysis of the process of 
halo formation (Vitvitska et al. 2002) shows strong evolution 
of the angular momentum of DM halos and its dependence 
upon the angular momenta of merged satellites. Perhaps, 
stronger correlation can be found between properties of the 
host cloud and the system of DM halos and its satellites. 

Let us note also that this scatter is partly enhanced by 
uncertainties in the determination of the complex shapes of 
the LSS elements. This is specially important for observa- 
tions when the shape of such elements can be found only 
approximately. 

5.6 The DM LSS elements and the Ly-a forest. 

The main observed characteristics of the Ly-a forest were 
analyzed in many publications (sec, e.g. Demianski et al. 
2006 and references there). Detailed discussion of the ob- 
served characteristics of the Ly-a absorbers can be found, 
for example, in Kim, Cristiani & D'Odorico (2002), Kim et 
al. (2002, 2004). 

There are some essential differences between the sample 
of simulated DM clouds and the observed Ly-a absorbers. 
The most important of them arc: 

(i) The Ly-a forest is related with the LSS elements with 
sizes, masses and other parameters much smaller than the 
simulated ones. 

(ii) The Ly-a forest is observed for larger neutral hydro- 
gen column density, Nhi > lO^^cm^^, which depends upon 
the UV background and is not directly linked with the DM 
component of absorbers. 

(iii) The evolution of the observed Ly-a absorbers de- 
pends upon the poorly known UV background and can sig- 
nificantly differ from the evolution of the simulated DM LSS. 

These differences prevent a direct quantitative compar- 
ison of characteristics of the simulated DM LSS and the ob- 
served Ly-a forest. However, the absorbers are dominated 
by the DM component, and, so, qualitative comparison of 
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their properties and evolution with the simulated one allows 
us to clarify some problems of the forest evolution. 

First of all the self similar character of evolution of the 
Doppler parameter, 6, is reflected in the surprising stabil- 
ity of its PDFs (Demianski ct al. 2006). The b parameter 
strongly depends upon the DM component and we can ex- 
pect that the most important aspects of its evolutions are 
closely correlated with the characteristics of evolution of DM 
clouds. 

The weak rcdshift dependence of the Doppler parameter 
and the high degree of matter relaxation along the line of 
sight for majority of the Ly-a absorbers are consistent with 
results obtained in Sec. 3.3, 3.5 and 4.1. It is important also, 
that the complex shape of the PDF P{b) is quite similar to 
the PDF P{a-s) plotted in Fig. 2. 

Our analysis in Sec. 4.2 cannot reproduce criteria of the 
Ly-a forest selection and in particular the redshift evolution 
of clumps separation (44, 45) does not represent the evo- 
lution of linear number density of the forest elements. The 
most important among the unknown factors are the disinte- 
gration of richer forest elements into a system of high den- 
sity subclouds linked by extended low density bridges, and 
variations of the ionizing UV background. Comparison of 
the LSS evolution in sirrmlations Siso and ^so confirms that 
the small scale perturbations accelerate the disintegration of 
the LSS elements and their internal structure becomes more 
complex. These problems need to be studied in more details 
with more representative and high resolution simulations. 
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